A new class of three-body states 
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We calculate the three-body spectrum for identical bosons interacting via attractive poten- 
tials. We have found an infinite number of three-body states even when the pair interactions are 
too weak to support any two-body states. These new states thus share this surprising scenario with 
the Efimov effect, but are not themselves Efimov states. Our effect occurs for both identical bosons 
and identical fermions, and it persists in the presence of two-body bound states. 
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In its long history, quantum mechanics has offered a 
wide variety of counterintuitive phenomena. One exam- 
ple that has received considerable interest lately is the 
Efimov effect from few-body physics. Predicted by Vi- 
taly Efimov in 1970 but not clearly observed exper- 

imentally until 2006 using ultracold atomic gases |5l4lO| . 
the effect that has now taken his name refers to the emer- 
gence of an infinite number of three-body bound states 
when none of the two-body subsystems are bound. Cen- 
tral to the effect is the requirement that the two-body in- 
teractions are short-ranged, i.e. fall off faster than 1 /r^ 
when r is a two-body interparticle distance, so that a 
two-body s-wave scattering length a can be defined. The 
Efimov effect occurs when |a|— >cx). 

Far from being an esoteric theoretical curiosity, the Efi- 
mov effect has a profound impact on low-energy three- 
body collisions even when a is finite 0-01 • lo'^ 
enough collision energies, the de Broglie wavelength is 
large compared to a characteristic short-range length 
scale To and thus insensitive to details on this scale. 
Consequently, in this limit, three-body systems display 
universal behavior dependent only on a [131] . Therefore, 
the desire to obtain the universal behavior of low-energy 
three-body observables motivates the study of weakly- 
bound three-body systems beyond the fundamental, in- 
trinsic interest such systems elicit. 

Generally speaking, two-body potentials can be cate- 
gorized as short- or long-ranged based on whether they 
fall off faster or slower, respectively, than 1/r^ when 
r—i'oo. Efimov's effect occurs in systems with short-range 
two-body interactions, and we know that an infinity of 
three-body bound states can occur in systems with long- 
range two-body interactions like the Coulomb interac- 
tion. Thus, to complete our knowledge of three-body sys- 
tems interacting via local two-body potentials, we must 
ask: What three-body spectrum results from attractive 
1/r^ two-body interactions — which lie exactly at the 
boundary between short- and long-range? 

This question turns out to be quite rich. For instance, 
depending on their strength, l/r^ potentials can have 
no two-body bound states when they are subcritical or 
critical, or they can have an infinity of two-body bound 
states when supercritical. Clearly, one interesting ques- 



tion is whether there is an infinite number of three-body 
bound states in the subcritical case. 

In this Letter, we show that this scenario is indeed 
possible for three interacting, indistinguishable bosons 
in three dimensions. We can thus reproduce the counter- 
intuitive result of Efimov's scenario, an infinity of three- 
body bound states in the absence of two-body bound 
states, but with long-range two-body potentials. While 
our scenario bears a superficial resemblance to Efimov's, 
it is, however, fundamentally different since a — a nec- 
essary ingredient for the Efimov effect — is not even de- 
fined. 

In addition, we show that when two-body bound states 
are formed, either by a short-range interaction or by 
strongly attractive 1/r^ interactions, an infinite number 
of three-body bound and resonant states appear. More- 
over, we show that our states also occur in a system of 
three identical fermions — a system for which no Efimov 
effect occurs. On this basis, we argue that the four-body 
"Efimov states" recently identified in Rcf [15|| are actually 
our states. 

Our study is based on numerical solutions of the 
three-body Schrodinger equation in the adiabatic hyper- 
spherical representation |16l4l8j . In this representation, 
the overall size of the system is measured by the hyperra- 
dius R, and the configuration of the system is represented 
by a set of five hyperangles collectively denoted il. Sub- 
stituting the exact wave function in this representation. 
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into the Schrodinger equation produces coupled hyperra- 
dial equations 



Wuu]f^+^Wuu'F^, = EF^,, (2) 



where /i = m/ \/Z for identical particles of mass to. The 
adiabatic channel functions ^i, in Eq. ([Ij are obtained 
from 
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in which the adiabatic Hamiltonian ffad inchides every- 
thing from the total Hamiltonian except the hyperradial 
kinetic energy. 

To determine the three-body spectrum, we focus on 
the behavior of W^i, in Eq. We will restrict our dis- 
cussion for three identical bosons to total orbital angu- 
lar momentum J—0 and total parity n=-|-l. Moreover, 
it is sufficient to neglect the interchannel coupling W^^' 
since the single-channel approximation gives a strict up- 
per bound to the three-body bound state energies [l9j. 
The are defined as 



W,,{R) = U,AR) + 



2n 



dR 



dR 
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Here, the double-bracket signifies that the integration is 
carried out only over the hyperangles. 

For the potential energy V that appears in ifad, 
we assume a pairwise sum of two-body potentials, 
V=v{ri2)+v{ri3)+v{r23), where 



1/4, 
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for interpartical distances . For the subcritical case, - 
l/4<a^<0 [iOl, v{rij) does not support any two-body 
bound states even though it is attractive and long- 
ranged. For the supercritical case, a^>Q, the singularity 
at the origin leads to the "fall-to-the-center" problem, 
and a ground state is not defined [il, After some 
regularization, v{rij) supports an infinite number of two- 
body bound states with properties that generically fol- 
low a geometrical scaling. For instance, the bound state 
energy En and mean radius (r)„ of the n-th state are, 
respectively, 



En+l/ En — e 



= „-27r/c 



and {r)n+i/{r)n = 
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For the two-body potentials in Eq. (I5|), it is well known 
that the hyperradial dependence can be exactly separated 
from the hyperangles. In this case, W^,^=Ui, and takes 
the form 



a' + 1/4 ^2 



2^iR^ 



(7) 



where are universal constants dependent only on a. 
For most any conceivable realization of an attractive 
two-body potential, however, there will be some 
short-range modification of Eq. ([5]) that will remove the 
singularity. To investigate the role of such a modification, 
we considered 
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1/4 



TO[rg sech^(rij/ro) + rfj] 



(8) 



The parameter tq sets the scale of the short-range regu- 
larization. 





FIG. 1: (Color online.) (a) The three-body effective poten- 
tials Ww(R) for three identical bosons with a^=Q. (b) The 
asymptotic behavior of WooiR) with the indicated values of 
a^. The symbols are solutions of Eq. ((3]), while the solid lines 
are fits based on Eq. 



Surprisingly, for the subcritical regime -1/4<q;^<0, 
the regularized potential, Eq. dramatically changes 
the asymptotic behavior of Woo from Eq. (O to 



Wqo{R) ^ - 



^l3ln{R/ro) + S ^ 



(9) 



Note, however, that unlike Eq. Q, Eq. © cannot cur- 
rently be justified rigorously. Instead, it is derived em- 
pirically from the numerical results. The quality of its 
asymptotic description is displayed in Fig. [U^b) where 
the axes have been manipulated so that Woo is a straight 
line if Eq. © holds. Moreover, by replacing sech."^ {rij /ro) 
in Eq. ([5]) by exp[— (rij/ro)^] and by unity, we have found 
that Woo from Eq. ^ is universal. That is, f3 is inde- 
pendent of the form of the short-range cutoff, but S is 
not. Equation © thus displays another departure of our 
effect from Efimov's since ro=0 and ro^O have the same 
asymptotic Wy,y(i?) in Efimov's case. 

While the R dependence of Woo in Eq. ^ is unusual, 
its consequence for the number of three-body bound 
states is clear: Woo falls off more slowly than l/R^, mak- 
ing it a long-range potential with an infinite number of 
three-body bound states. Figure [ija), which shows the 
lowest several W,y^ for the most attractive case a^=0, 
also makes it clear that trimers produced by this effect 
are going to be very weakly bound and very large — the 
minimum in Woo lies at roughly i?=7000ro. As de- 
creases, Woo becomes shallower, and the position of the 
minimum moves towards even larger R. 

Unfortunately, the R dependence of Woo in Eq. ([9|) 
does eliminate the possibility of finding an exact result 
like Eq. © for the three-body spectrum. But, given the 
extremely slow variation of the numerator in Eq. ^ and 
the fact that the hyperradial probability density is con- 
centrated at the outer turning point, we can treat Woo as 
a pure potential with a variable coefficient based on 
the size of each state to find an approximate expression 
for En- Replacing R in the numerator of Eq. ^ by 
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FIG. 2: (Color online.) (a) Comparison of the three-boson 
bound state energies calculated from Eq. ([2]) (circles) and 
from Eq. (|10|) (triangles). The lines are added to guide the 
eye. (b) Fits of the potential parameters /3 and 5 from Eq. ((9|: 
symbols denote individual fits of Eq. (|9]) to the numerical Woo; 
and solid lines, fits of Eq. pip to these values. 



and using Eq. (jS]), we find for n3>l 




FIG. 3: (Color online.) (a) The lowest four three-body ener- 
gies calculated from Eq. (I2J for Woo from Eq. with a hard 
wall added at i?=100ro. The dashed line denotes the lowest 
two-body threshold, (b) Schematic three-body binding ener- 
gies assuming there is a deeply-bound two-body state whose 
energy is independent of q. 



En+i/En = exp - 
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(10) 



having neglected 5. 

Figure [2][a) compares the three-body spectrum given 
by Eq. (ITU|) with the energies calculated by solving Eq. ^ 
numerically using the potential given by Eq. © cut off 
by a hard wall at R—\QQro. The agreement between the 
two is quite good, especially given that the energy scale 
covers 60 orders of magnitude, underscoring the value of 
Eq. ()10|) since calculating many these states numerically 
is actually quite challenging. The /3 used in Fig. 
were determined by fitting the empirical expressions 

/3 = exp(6^a^) -f c,3, —5 — aseyi\){bsa'^) + cs- (11) 

to the numerical data as shown in Fig. [UJb). 

An infinity of three-body bound states thus exists at 
a^=0, and Fig. EJb) shows that an infinity also exists 
for some negative . Figure [2l[b) further shows that /3 
crosses zero at a critical value aj. that can be determined 
from Eq. ^ to be a^«-0.0072. However, due to the 
difficulty in solving Eq. ([3]) numerically to the large dis- 
tances required to identify the asymptotic behavior — at 
least i?=10^ro near a1 — this critical value should be 
treated as an approximate upper bound. The true value 
of may be even more negative. We have thus found 
an infinite number of three-body states for a range of 
where no two-body states exist. 

We summarize the behavior of the three-body spec- 
trum discussed so far in the a^<0 portion of Fig.[3ja). At 
a^=aj., an infinite number of three-body states appears 
with energies given by Eq. (|10p . At a^=0, an infinite 
number of two-body states emerges — only the lowest of 
which is indicated in Fig. [Sja) for simplicity. 

Before considering the a^>0 behavior, let us revisit 
the zero-range two-body potential in Eq. The theo- 
retically attractive fact that R can be exactly separated 



in the Schrodinger equation begs the question: what are 
the al in Eq. dZ])? It turns out that Eq. © at a fixed 
R provides the answer. In the limit rg— >0, the numer- 
ator — proportional to — diverges. This conclusion 
is supported by taking the limit numerically for various 
regularization schemes. While this divergence is a bit 
disappointing, it is also interesting in its own right as 
a novel three-body "fall-to-the-center" problem that oc- 
curs even when there is no divergence in the two-body 
problem. 

Returning to Fig. |31[a), we see that the three-body en- 
ergies connect smoothly across a^=0. The nature of the 
spectrum, however, changes. For a^>0, the three-body 
system can break up into a dimer and a free particle 
when i?— >-cx). The leading term in W^i, is then given by 
the sum of the two-body interactions between the free 
particle and each of the particles in the dimer. Thus, the 
asymptotic form for the J^={)'^ three-body potentials is 



vl 



(12) 



where Eyi is the energy of a two-body vibrational state 
V with orbital angular momentum I and 



—a 
3 



;(/ + 1). 

12 ^ ' 



(13) 



Note that to have J'^^O^, / must also be the orbital 
angular momentum of the free particle relative to the 
dimer. We have verified Eq. (jl2p numerically for the 
two-body potential in Eq. ([5]). 

For s-wave dimers, Eq. (|13p shows that for each two- 
body bound channel, the interaction in Eq. (|12p is al- 
ways supercritical when a^>0. An infinite number of 
three-body bound states thus lie in the hyperspherical 
potential asymptoting to Eqo with energies relative to 
Eqq given by Eq. © after replacing a by acs- But, there 
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are also an infinite number of three-body resonances at- 
tached to each of the excited two-body thresholds E^o 
whose energies relative to their respective threshold are 
also given approximately by Eq. 

For larger values of a^, it is possible to have an infinity 
of two-body dimers with Z>0. In particular, when /o('o+ 
l)<a^<{lo + l){lo + 2), the effective two-body potential, 
including the centrifugal barrier, is supercritical for l<lo. 
Equation ([T3| shows that so long as is supercritical for 
a given I, a^g will also be supercritical. In other words, 
there will be infinite series of states attached to every 
two-body threshold. While the energies of all of these 
series relative to E^^i are given by Eq. (jS]), note that acff 
in Eq. ([T^ does depend on I. 

To this point, we have assumed that the short-range 
modification of the two-body potential only makes it 
shallower, thus preserving the property that no two- 
body states exist for a^<0. This assumption is conve- 
nient for establishing the characteristics of the system 
but is certainly not necessary. If we now allow short- 
range modifications that support deeply-bound two-body 
states, then by the arguments above, it is clear that 
Eq. (IT2|) will also apply for < 0. Because these deeply- 
bound two-body states are independent of a, Eq. im- 
plies a new critical value for determined from 
ajj=3l{l+l)/8—5/32. The spectrum of these three-body 
states is given by Eq. and they will coexist with 
those states from the higher channel given by Eq. ([T0|. 
We show the lowest few states of the three-body spec- 
trum attached to this deeply bound two-body threshold 
in Fig. ^h). There will be a similar spectrum for each 
deeply bound two-body state. 

A natural question to ask is whether the effect we have 
described occurs for other systems such as three identical 
bosons with J>0, different mass particles, or even com- 
binations of identical fermions. Since the latter seems to 
have some relevance to the work recently published in 
Ref. [31, we will focus on it. 

We have calculated W„i, numerically for three indistin- 
guishable, spin-polarized fermions interacting via Eq. ([8|) . 
In this case, we study the symmetry J^~l^ since it gives 
the lowest potential in the non-interacting case 22] ■ The 
lowest two-body angular momentum satisfying the sym- 
metry constraints is 1=1, implying that can be as large 
as 2 without having a two-body bound state. Having ex- 
amined Woo for several a^<2, we have empirically found 
that the asymptotic behavior of Woo is 



WooiR) 



2^1n(i?/ro)i?2 



(14) 



For a^=2, we have obtained a^fj =5.24 and 7=4.19 by fit- 
ting Eq. (|14p to our numerical potential. Therefore, our 
effect does occur: the three-fermion system has an infi- 
nite number of three-body bound states — in the absence 
of two-body bound states — with a spectrum given ap- 
proximately by Eq. ^ . Note that even though our effect 



shares the same spectrum as the Efimov effect, it exists 
in a 1"*" system of identical fermions where the Efimov 
effect does not [23| . 

Based on these results, we believe that our effect might 
better explain the four-body states recently found by 
Castin et al. 15| than does the Efimov effect. Castin 



et al. considered a system of three heavy, spin-polarized, 
identical fermions H plus a fourth lighter particle L, as- 
suming the only interactions were contact potentials be- 
tween H and L. Upon setting the H+L s-wave scattering 
length ghl to infinity, they found a supercritical in- 
teraction for 13.384<m/f/TOL<13.607 — if the Hs had 
1"*" symmetry. 

To connect with our effect, we can approximately re- 
duce their four-body problem to an effective H+H+H 
problem by integrating out the L motion via the Born- 
Oppenheimer approximation (j^ ]. The resulting Born- 
Oppenheimer surface for most configurations of the three 
Hs can be well-approximated as a sum of H+H pair po- 
tentials. Since a//L=oo, Efimov's analysis tells us these 
pair potentials behave as l/r^ [l3| . We thus have a three- 
body problem with attractive l/r^ interactions — i.e. 
the focus of this Letter. 

Equation (fT4|l shows that, unlike the bosonic case, we 
can take the limit rp— »0 for fermions to match the as- 
sumptions in Ref. [la |. leaving just the first term in 
Eq. (fTi|) in agreement with Eq. ([7]). For a^=2, corre- 
sponding to m^f/mL=13.607, we thus obtain for aQ from 
Eq. ([7]) the supercritical value of 5.24. Reducing un- 
til ckQ is zero gives a^^l.G — the fermion equivalent of 
found for bosons above. Remarkably, we have thus 
found an infinity of three-body bound states when the 
effective two-body potential (including the centrifugal 
term) is repulsive! Using the zero-range model [ll,[l3, 
we can convert this to a mass ratio of mn /'mL=11.58, 
implying that an infinity of 1+ three-fermion bound 
states — in the absence of two-body states — exists for 
11.58<m///mL<13.607. 

It is clear that our range of mass ratios and values of 
a^g do not quantitatively match those from Ref. [l5|. 
One of reasons for these quantitative differences is likely 
the fact that m/f/rTiL is only roughly 10, reducing the 
reliability of the Born-Oppenheimer approximation. An- 
other likely reason for the difference is the approximation 
of the Born-Oppenheimer surface as purely a pair-wise 
sum of 1 /r^ interactions when simple arguments suggest 
that the surface is less attractive for some configurations. 
Nevertheless, we think our explanation of the infinity of 
states from Ref. is compelling. And, while they could 
not exist without the Efimov physics of the H2L system, 
they are not themselves Efimov states for all of the rea- 
sons we have discussed. 

In summary, we have found a new class of three-body 
states: when three particles interact via attractive l/r^ 
potentials, they can have an infinite number of three- 
body bound states even when there are no two-body 
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bound states. This effect occurs for both 0+ identical 
bosons and 1+ identical fermions — and almost certainly 
for many other cases as well. For identical bosons, this 
effect produces a unique three-body spectrum and a new 
kind of three-body fall-to-the-center problem. For identi- 
cal fermions, this effect sounds especially surprising since 
it occurs when the effective two-body interaction is com- 
pletely repulsive. Moreover, this effect still occurs when 
there are two-body bound states. While these states do 
share some characteristics with Efimov states, their phys- 
ical origin is quite distinct. We use this distinctiveness to 
argue that the recently-found four-body "Efimov" states 
are, in fact, "our" states — arguments that seem likely 
to apply equally well to even more particles. 
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